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Abstract. In this work, we study the uncertainty of nuclear model parameters for 
neutron induced 56 Fe reactions in fast neutron region by using the Total Monte Carlo 
method. We perform a large number of TALYS runs and compare the calculated 
results with the experimental data of the cross sections to obtain the uncertainties of 
the model parameters. Based on the derived uncertainties another 1000 TALYS runs 
have been performed to create random cross section files. For comparison with the 
experimental data we calculate a weighted \ 2 value for each random file as well as the 
ENDF/B-VII.l, JEFF3.1, JENDL4.0 and CENDL3.1 data libraries. Furthermore, we 
investigate the optical model parameters correlation obtained by way of this procedure. 
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1. Introduction 

Since the beginning of the century, the nuclear science community is putting more 
and more attention to the assessment of uncertainties which is important for both basic 
physics and technological applications. More recently, a new method has been developed 
and applied based on Monte Carlo calculations and called "Total Monte Carlo", or 
TMC[lJ. This method relies on a large number of nuclear model calculations, with 
different random nuclear data in each of them. The essential idea of this method is to 
assume that each nuclear model parameter has its own uncertainty. Model calculation 
will be performed many times, whereby each time all elements of the input parameter 
vector are randomly sampled from a normal or uniform distribution, centered around a 
initial value, with a specific width for each parameter. At NRG, this method has been 
implemented in a code called TASMAN[2J. In TASMAN, a binary reject/accept method 
is used to introduce correlations between the nuclear model parameters: calculations 
are judged by their agreement with experimental data(see sect. II A) and a TALYS [3] 
run with a certain parameter set might be either accepted or rejected. In this way 
only certain model parameter combinations survive, and the parameter distribution is 
automatically determined numerically by experimental data, without having to resort 
to a priori distribution. 

By using this method we investigate the uncertainty of model parameter and 
parameter correlations for n+ 56 Fe reactions. Based on the derived parameter 
uncertainties another a set of 1000 random cross section files has been created with 
TALYS. For comparison with the experimental data we calculate a weighted \ 2 value 
for each random file as well as the existing evaluated data libraries, such as ENDF/B- 
VII.l, JEFF3.1, JENDL4.0 and CENDL3.1. Section II describes the method in more 
detail and gives the main results. A briefly summary is given in last section. 

2. Methodology and results 

2.1. Methodology 

We use the following scheme in our work. 

• Using the TALYS default parameters as starting point a manual search for a good 
set of model parameters that well describe experimental data in the considered 
energy region is performed. Note that resonances are not included in this work. All 
experimental data used in this work are taken directly from the EXFOR database. 

• For the next step we make a rough assumption on parameter uncertainties. These 
are used to sample parameter space in a number of TALYS runs. Given typical 
experimental uncertainties a TALYS run is accepted if it does not differ from 
the average result by more than the experimental uncertainty, and it is rejected 
otherwise. In this way we find both improved values for parameter uncertainties 
(given in Table [l| and obtain correlations between model parameters (see Sect. 
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Table 1. Uncertainty of some nuclear model parameters for 56 Fe, given by fraction(%) 
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II.C). 

• Finally, using the derived parameter uncertainties, a set of 1000 randomized 
cross section files is obtained from TALYS and this set is directly compared to 
experimental data by means of a x 2 test as described below. 

The run with the lowest x 2 value together with the parameter uncertainties derived 
in the fashion described above are a good starting point for use in TMC. 

2.2. x 2 calculation 

we calculate x 2 values for comparison with experimental data. First, the average xl f° r 
each reaction channel is calculated as follows 

1 N a 1 - a 1 2 



Where N is the number of experimental data of each reaction channel, gt the calculated 
result, ge the experimental data, Age the uncertainty of the experimental data. The 
total average x 2 f° r eacn random file is defined as a sum over all reaction channels N c 

X = • (2) 

Ei Wi 
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Table 2. List of 56 Fe experimental data selected for x 2 calculation with their weight. 



Reaction channel 


Energy points 


Energy region ( Me V) 


weight 
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(n,2n) 
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11.88-19.3 


11.6 
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We weight the x 2 with 



Wi = y ~^=Na max . (3) 

Here Acfg is the average uncertainty of each reaction channel, a max the maximum value 
of the cross section. A large weight has been given for the reaction channel with a large 
number of experimental data, more precise data and large cross section. By this way, 
the total, elastic, nonelastic, inelastic reaction with large number data and large cross 
section, play the most important roles in \ 2 calculation (see Table [2|. 

Fig [I] presents the results of 1000 random 56 Fe files in terms of x 2 as defined in 
Eq|2j where each random x 2 is represented by an entry into the histograms. We also 
compare to the existing evaluations including ENDF/B-VII.l, JEFF3.1, JENDL4.0 and 
CENDL3.1. It can be seen from the figure that a number of our random cross section 
files perform as well as and in some cases even better than the evaluated libraries. As 
can be seen, the distribution is not symmetric, and has a tail towards high x 2 values jl]. 



2.3. Parameter correlations 

The parameter correlations have been investigated by using the TASMAN code. In this 
case, we select any two of optical model parameters to sample with large uncertainty 
ranges while keeping the other parameters constant. Fig [2] shows the accepted and 
rejected cases decided by r v against v\. We use the larger uncertainties in this study in 
order to show the boundary of the accepted region. Considering the relation between 
r v and v%, which is v\ *r% = constant, r v and v\ show the expected negative correlation. 
Actually, most optical model parameters show similar behaviour for the cross sections 
such gainst a v and d\ against aj, etc, which have negative correlations. Fig 
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Figure 1. \ 2 values distribution for the random 56 Fe files per bin, compared to the \ 2 
values for the ENDF/B-VII.l, JEFF3.1, JENDL4.0 and CENDL3.1 libraries marked 
by the vertical lines with the corresponding color showed in figure. 



[3] shows corresponding results for a v against t> 2 , and indicates that a v and v 2 have a 
positive correlation, i.e. a v and v 2 have opposite effect on the cross section. It can be 
seen from the figures that the slope of the accepted area is a measure for the parameter 
correlations. 

3. Summary 

In this work we obtain a set of uncertainties for the most important model parameters 
for n+ 56 Fe reactions. 1000 TALYS runs have been performed to create random cross 
section files based on the derived uncertainties. For comparison with the experimen- 
tal data, we calculate \ 2 values with different weights for different reaction channels. 
Larger weights have been given for the reactions with a large number of experimen- 
tal data, more precise data and large cross section. Judging from the obtained x 2 > a 
number of random files perform as well as the evaluated libraries. Parameter corre- 
lation have been investigated by using TASMAN code. The results show correlations 
between optical model parameters. Meanwhile, this study further verified the reliability 
of the TASMAN code and the described reject/accept method. Actually, the whole 
TMC method include reactor code simulations, i.e. bechmark testing process. We also 
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Figure 2. Accepted and rejected cases decided by r v against v\ 

use TMC method for reactor code simulations on burn-up [5 J and safety parameters [6J. 
Next step we will produce the ENDF random files including resonance part by using 
the whole TASMAN package for reactor code simulation. 
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Figure 3. Accepted and rejected cases decided by a v against V2- 



